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We review the contribution of extra galactic inspiralling double neutron stars, to the LISA astro- 
physical gravitational wave foreground. 

Using recent fits of the star formation rate, we show that sources beyond z, — 0.005 contribute 
to a truly continuous background, which may dominate the LISA instrumental noise in the range 
3 ~ 10 -4 - 1 x 10~ 2 Hz and overwhelm the galactic WD- WD confusion noise at frequencies larger 
than c ~2x 10~ 3 . 
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& • I. INTRODUCTION 

Q 

Compact neutron star binaries are among the most promising sources of gravitational waves. At low frequencies, 
the continuous inspiral signal may be detectable by the space antenna LISA, while ground based interferometers such 
as VIRGO 0, LIGO @, GEO Q or TAMA ff|, are expected to detect the last few minutes prior coalescence, at 
frequencies up to 1.4-1.6 kHz. 

In a first paper Q (hereafter paper I), we have investigated the high frequency signal and its detection with the 
first generations of ground based detectors. Using new estimates of the mean merging rate in the local universe, 
that account for the galactic star formation history derived directly from observations and include the contribution of 
t-H ' elliptical galaxies, we predict a detection every 148 and 125 years in the volume probed by initial VIRGO and LIGO, 
, and up to 6 detections per year in their advanced configuration. 

In a second paper 0] (hereafter paper II), we used numerical simulations to estimate the gravitational wave emission 
produced by the superposition of unresolved extra-galactic sources. As in paper I, we were interested in the few sixteen 
^ . minutes before the last stable orbit is reached, when more than 96% of the gravitational wave energy is released and 
Qh when the frequency evolves in the range 10-1500 Hz, covered by ground based interferometers. We find that above the 
critical redshift z* = 0.23, the sources produce a truly continuous background, with a maximal gravitational density 
parameter of £l gw ~ 1.1 x 10~ 9 around 670 Hz. The signal is below the sensitivity of the first generation of ground 
based detectors but could be detectable by the next generation. Correlating two coincident advanced-LIGO detectors 
C$ • or two third generation EGO interferometers, we derived a S/N ratio of 0.5 and 10 respectively. 

In this article, we extend our previous simulations to investigate the stochastic background from the early low 
frequency inspiral phase. Estimates of the emission from the various population of compact binaries (see for instance 
[cl M, [Tol. fill. H3, [I3T| for the galactic contribution and [III, [ijl [l6| for the extra-galactic contribution) or from captures by 
super massive black holes ([I3])> which represent the main source of confusion noise for LISA, are of crucial interest for 
the development of data analysis strategies. The paper will be organized as follow. In II., we describe the simulations 
of the DNS population, in III. the contribution of inspiralling DNS to the stochastic background is calculated and 
compared to the LISA instrumental noise and to the WD- WD galactic foreground. Finally, in IV. the main conclusions 
are summarized. 



II. SIMULATIONS OF THE DNS POPULATION 



To generate a population of DNSs (each one characterized by its redshift of formation and coalescence timescale) 
we follow the same procedure as in paper II, to which the reader is referred for a more detailed description. 
The redshift of formation is randomly selected from the probability distribution [l8j : 

P/(*/) = ^ (1) 
Hp 
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constructed by normalizing the differential DNS formation rate, 

«*(.,) M 

The normalization factor in the denominator corresponds to the rate at which massive binaries are formed in the 
considered redshift interval, e.g., 

R p = [ {dRf(zf)/dzf)dz f , (3) 
Jo 

which numerically gives R p = 0.044 s" 1 in the two cosmic star formation rate adopted in this paper. 
The formation rate of massive binaries per redshift interval eq. [5] writes : 

R*^f)-\ 1 + Zf dZf (4) 

where R*j(z) is the cosmic star formation rate (SFR) expressed in M Mpc _3 yr _1 and X p is the mass fraction converted 
into DNS progenitors. Hereafter, rates per comoving volume will always be indicated by the superscript '*', while 
rates with indexes "z/" or "z c " refer to differential rates per redshift interval, including all cosmological factors. The 
(1+z) term in the denominator of eq. [4] corrects the star formation rate by time dilatation due to the cosmic expansion. 
The element of comoving volume is given by 

with 

E(a l ,z) = [n m (i + zf + a v } 1 / 2 (6) 

where f2 m and £l v are respectively the present values of the density parameters due to matter (baryonic and non- 
baryonic) and vacuum, corresponding to a non-zero cosmological constant. Throughout this paper, the 737 flat 
cosmology [19], with fl m = 0.30, Q v — 0.70 and Hubble parameter Ho = 70kms~ 1 Mpc~ 1 c 2Q] is assumed. 

Recent measurements of galaxy luminosity function in the UV (SDSS, GALEX, COMB017) and FIR wavelengths 
(Spitzer Space Telescope), after dust corrections and normalization, allowed to refine the previous models of star 
formation history, up to redshift z ~ 6, with tight constraints at redshifts z < 1. In our computations, we consider 
the recent parametric fits of the form of [2l| , provided by [22| , constrained by the Super Kamiokande limit on the 
electron antineutrino flux from past core-collapse supernovas. It is worth mentioning that the final results are not 
significantly different if we adopt, as in paper II, the SFR given in (23|. 

Throughout this paper, we assume that the parameter A p does not change significantly with the redshift and can 
be considered as a constant. In fact, this term is the product of three other parameters, namely, 

Ap = ^Nsfb^NS (?) 

where /3jvs is the fraction of binaries which remains bounded after the second supernova event, fb is the fraction of 
massive binaries formed among all stars and Xns is the mass fraction of neutron star progenitors. From paper I we 
take_B NS = 0.024 and f b = 0.136. 

[22j investigated the effect of the initial mass function (IMF) assumption on the normalization of the SFR. They 
showed that top heavy IMFs are preferred to the traditional Salpeter IMF [24[ and their fits are optimized for IMFs 
of the form: 

„ , t^" 1 - 5 for 0.1 < to < m 
v ' ^ for m < to < 100 

with a turnover below too = 1 M Q , and normalized within the mass interval 0.1 - 100 M such as J m£(m)dm = 1 



In this paper, we assume 7 = 2.35 (A modified Salpeter), but taking 7 = 2.2 ( 25|) would not change the final 
results significantly. To be consistent with the adopted models of the SFR, we follow 22j and assume a minimal initial 
mass of 8 M Q for NS progenitors, the upper mass being 40 M©. It results finally , Xns = J 8 £(m)dm — 9 x 10~ 3 
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The next step consists to estimate the redshift z b at which the progenitors have already evolved and the system is 
now constituted by two neutron stars. This moment fixes also the beginning of the inspiral phase. If r b (« 10 8 yr) is 
the mean lifetime of the progenitors (average weighted by the IMF in the interval 8-40 M©) then 

z b = z f - H r b (l + z f )E(z f ) (9) 

Once the beginning of the inspiral phase is established, the redshift at which the coalescence occurs should now 
be estimated. The duration of the inspiral depends on the orbital parameters and the neutron star masses. The 
probability for a given DNS system to coalesce in a timescale r was initially derived by [261 ] and confirmed by 
subsequent simulations 0, [23] and is given by 

P(t) = B/t (10) 

Simulations indicate a minimum coalescence timescale tq = 2 x 10 5 yr and a considerable number of systems having 
coalescence timescales higher than the Hubble time. The normalized probability in the range 2 x 10 5 yr up to 20 Gyr 
implies B = 0.087. Therefore, the redshift z c at which the coalescence occurs is derived by solving the equation 

r = r(z c ,z b ) (11) 

where 

dz 



z)E(z) 



(12) 



III. THE GRAVITATIONAL WAVE BACKGROUND 

Compared with the previous study, slight complications arise when calculating the spectral properties of the signal. 
In paper II, we were interested in the signal emitted between 10-1500 Hz, which duration in the source frame is short 
enough, less than 1000 s, to be considered as a burst located at the redshift of coalescence. Now, as we are looking at 
inspiral phase in the LISA frequency window (10~ 4 — 0.1 Hz), the evolution of the redshift of emission with frequency 
must be taken into account. For the same reason, the critical redshift at which the population of DNSs constitute a 
truly continuous background can't be determined by simply solving the condition D(z) > 1, where D is the duty cycle 
defined as the ratio of the typical duration of a single burst f to the average time interval between successive events 
(see paper II). Instead, we remove the brightest sources, lying in the close Universe, where the distribution of galaxies 
is expected to be highly anisotropic, in particular due to the strong concentration of galaxies in the direction of the 
Virgo cluster and the Great Attractor (see paper I), and where the density is well below the average value derived 
from the SFR (f29[). In a second step, we verify that the resulting background is continuous, by using both Jaque 
Bera and Lilliefors gaussianity tests with significant level a — 0.05 over a set of 100 realizations. 

For each realization, the number of simulated DNSs corresponds to the expected number of extragalactic DNSs 
observed today and is derived as, 

N p = [ N(z)dz (13) 
Jo 



z' b and z. Thus 



The DNSs present at redshift z were formed at z' b > z, with a coalescence time larger than the cosmic time between 

,6 

N(z) = \ R zc {z')r,{z')dz' (14) 



where 

r)(z')= i P(r)dT (l.r,. 



20 Gyr and T mi „=Max(2 x 10 5 yr ;r(z, z' b )), which gives N p ~ 1.4 x 10 6 . 
We obtain that after z* = 0.005, corresponding to the distance beyond the Virgo cluster (~ 25 Mpc), when the 
relative density fluctuations are < 0.5 ([29J]), the background is continuous in the range 10~ 4 — 0.1 Hz. We consider also 
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a more conservative threshold of z* = 0.02, corresponding to the distance beyond the Great Attractor (~ 100 Mpc), 
when the relative density fluctuations are < 0.05 and when the Universe is expected to become homogeneous ( [29i |30T] ) . 
The gravitational flux in the observer frame produced by a given DNS is: 

where z e is the redshift of emission, cLl = (1 + z e )r the distance luminosity, r the proper distance, which depends 
on the adopted cosmology, dE gw /dv the gravitational spectral energy and v = (f + z e )v the frequency of emission in 
the source frame. 

In the quadripolar approximation and for a binary system with masses mi and 11J2 in a circular orbit: 

dE gw jdv = Kv- 1/3 (17) 
where the fact that the gravitational wave frequency is twice the orbital frequency was taken into account. Then 

(Gtt) 2 / 3 mi m 2 , , 

K = q 7 1 u73 

3 (mi+m 2 ) L/ - i 

Taking mi = m 2 = 1.4 M Q , one obtains K = 5.2 x 10 50 erg Hz -2 / 3 . 
The redshift of emission is derived by solving the relation [l4| : 

is- 8/3 = is L so + K'{T(z b ,z e )) (19) 

where 

K' = <^!^!1 m 

5c 5 v ; 

(or numerically K' = 67.83 yr _1 ), and where the time spent between redshifts z e and z;, is given by: 

r° dz . , 

r(z b ,z e =t- / — (21) 

A. Properties of the background 

The spectral properties of the stochastic background are characterized by the dimensionless parameter [28[ : 

ttgw{Vo) = -^v F Vo (22) 

where v Q is the wave frequency in the observer frame and p c the critical mass density needed to close the Universe, 
related to the Hubble parameter H by, 

* = 8^ < 23 ) 

F Vo is the gravitational wave flux at the observer frequency (given here in erg cm _2 Hz _1 s~ 1 ) v 0) integrated over 
all sources above the critical redshift z*, namely 



F Uo = / fuJRc(z) (24) 

In our simulations, the integrated gravitational flux is calculated by summing individual fluences feq. I16[) . scaled 
by the ratio between the total formation rate of progenitors (eq. [3]) and the number of simulated massive binaries, as: 

n p 

y i—l 

Figure [1] shows the density parameter VL gw as a function of the observed frequency averaged over 10 simulations. The 
continuous line corresponds to sources located beyond the critical redshift to have a continuous background, z* = 0.005 
(or d ~ 25 Mpc), while the dashed line corresponds to the more conservative estimate, with a cutoff at z* = 0.02 (or 
d ~ 100 Mpc). After a fast increase below ~ 5 x 10~ 5 Hz, Cl gw increases as v^ 3 to reach fl gw ~ 3 x 10 -10 at 0.1 Hz 
for z* = 0.005 (lip,,, ~ 10~ 10 for z* = 0.02). This is one order of magnitude larger than the previous predictions by 
Schneider et al. [14| . 
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FIG. 1: Spectrum of the expected gravitational energy density parameter Q gm , corresponding to NS-NS inspirals occurring 
beyond d ~ 25 Mpc (continuous curve) and d ~ 100 Mpc (dashed curve). 



B. consequences for LISA 

Astrophysical backgrounds represent a confusion noise for LISA, which spectral density is given by [T^: 

sr fe K) = S^ g wK) (26) 

Figure [2] shows S^ ck {v ) compared to the LISA sensitivity [33| and to the contribution of unresolved galactic WD- WD 
P |33| . The contribution from DNSs may dominate the LISA instrumental noise between ~ 3 x 10~ 4 - 10~ 2 Hz for 
= 0.005 (~7x 10~ 4 - 6 x 10~ 3 Hz for z* = 0.02) and the galactic WD- WD confusion noise after ~ 2 x 10~ 3 . 
However, the resulting reduction in the sensitivity should be less than a factor 4 and thus shouldn't affect significantly 
signal detection. 

The critical redshift at which a single DNS can be resolved and detected by LISA (Figure [3]) is at least one order 
of magnitude smaller (z m i n ~ 3 x 10~ 4 at frequencies v a > 0.01) than our threshold to have a gaussian background 
= 0.005 (or z* = 0.02) . Sources lying between the two regimes are responsible for a non gaussian and non isotropic 



cosmic "popcorn" noise (see 3J] for a recent review) that will be investigated in a further work. This contribution, 
which statistical properties differ from those of both the instrumental noise and the cosmological background, could be 
detected and removed by new data analysis techniques currently under investigation, such as the search for anisotropies 
HH that can be used to create a map of the GW b ackg round [3(|[37|, the maximum likelihood statistic [Hj], or methods 
based on the Probability Event Horizon concept [39| , which describes the evolution, as a function of the observation 
time, of the cumulated signal throughout the Universe. 
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FIG. 2: Gravitational strain in Hz corresponding to NS-NS inspirals occurring beyond d ~ 25 Mpc (continuous curve) 
and d ~ 100 Mpc (dashed curve), along with the LISA instrumental noise (black) and the galactic WD- WD foreground (light 
gray). 



IV. CONCLUSIONS 



In this paper, we have modelled the contribution of extra galactic inspiralling double neutron stars, to the to the 
astrophysical gravitational wave foreground, in the range of sensitivity of the space detector LISA. 

Using a recent fit of the star formation rate [22], optimized for a top heavy A modified Salpeter IMF, we show that 
sources beyond z* = 0.005 (or z* = 0.02 for a more conservative estimate) constitute a truly continuous background. 
We find that the density parameter Q gw , after a fast increase below ~ 5 x 10 -5 Hz, increases as v^ 3 to reach 



il gw ~ 3 x 10 -10 (Ct gw ~ 10 -10 ) at 0.1 Hz, which is one order of magnitude above the previous estimate of [14[. As 
a result, the signal may dominate the LISA instrumental noise in the range ~ 3 x 10~ 4 - 10 -2 Hz (~ 7 x 10~ 4 - 
6 x 10 -3 ) and overwhelm the galactic WD- WD confusion noise at frequencies larger than v a ~ 2 x 10~ 3 . Sources 
located closer than z*, but still too far to be detected by LISA, constitute a non gaussian and non isotropic cosmic 
"popcorn" noise, which will be investigated in a further work. 
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FIG. 3: Maximal redshift as a function of the observed frequency, at which a single DNS can be detected by LISA after one 
year of integration and with a signal-to-noise ratio S/R = 5. The galactic WD- WD foreground is taken into account in the 
LISA sensitivity. 
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